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Математическое моделирование процессов и явлений часто заменяет с высокой степенью достоверности дорогостоящие экспе- 
риментальные исследования, особенно в такой области как физика плазмы и плазмохимия. Экспериментальные исследования 
физико-химических процессов, протекающих в многокомпонентных плазменных потоках, например, при взаимодействии воз- 
душной низкотемпературной плазмы с теплозащитными материалами и пиротехническими составами, требуют предваритель- 
ной информации о характеристиках таких плазменных систем. Например, при масс-спектрометрических исследованиях соста- 
ва многокомпонентной плазмы важно предварительно знать, какие компоненты могут образовываться в результате протекаю- 
щих в плазме процессов. С этой целью для оценки состава многокомпонентной плазмы были использованы методы термодина- 
мического анализа. 


Известно, что в пиротехнических изделиях, 
применяемых с целью излучения в инфракрасной 
области спектра, используются в качестве окисли- 
телей нитрат натрия в сочетании с нитратом цезия 
(или нитратом калия), а в качестве горючего - по- 
рошок магния. Каждый из компонентов пиротех- 
нических составов (ПТС) вносит определенный 
вклад в протекающие процессы при их взаимодей- 
ствии с воздушной плазмой. 

В связи с этим возникает необходимость иссле- 
дования процессов взаимодействия индивидуаль- 
ных составляющих пиротехнических составов с 
плазмой, что позволяет с большей достоверностью 
оценивать механизмы процессов взаимодействия 
плазмы с ПТС. Кроме того, знание механизмов 
процессов взаимодействия плазмы с нитратами 
щелочных металлов (НЩМ) позволяет повысить 
эффективность пиротехнических изделий. 

Учитывая многообразие компонентов, образу- 
ющихся в результате физико-химических процес- 
сов при взаимодействии плазмы с веществом, не- 
обходима предварительная информация о возмож- 
ном выходе продуктов реакций и зависимости их 
концентрации от температуры и давления в плазме. 
С учетом этого, перед проведением эксперимен- 
тальных исследований проводился термодинами- 
ческий анализ исследуемых сложных многокомпо- 
нентных плазменных систем. 


Метод расчета основан на использовании прин- 
ципа максимума энтропии изолированной термо- 
динамической системы. Это связано с тем, что на- 
иболее общим параметром, по которому можно су- 
дить, находится ли термодинамическая система в 
состоянии равновесия или нет, является величина 
энтропии системы, поскольку энтропия связана со 
степенью упорядоченности энергетического состо- 
яния микрочастиц, из которых состоит рабочее те- 
ло (или исследуемый объект). 

Данный метод позволяет проводить расчет рав- 
новесного состояния сложных многокомпонент- 
ных систем, включающих не только индивидуаль- 
ные газообразные химические вещества в электро- 
нейтральном и ионизированном состояниях, но и 
конденсированные вещества и их растворы. 

Расчет параметров равновесного состояния 
сложных многокомпонентных плазменных сред по 
принципу максимума энтропии проводится по сле- 
дующей системе уравнений [2]: 
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Известно, что определение равновесного соста- 
ва сводится в большинстве методов к совместному 
решению системы уравнений закона действующих 
масс, уравнений материального баланса и уравне- 
ния нормировки (уравнения Дальтона). 

Однако большинство созданных программ для 
решения системы уравнений носит частный харак- 
тер, что ограничивает их применение для расчета 
равновесия при произвольном наборе химических 
элементов. 

В связи с этим рассмотрим наиболее совершен- 
ный метод расчета равновесного состояния слож- 
ных многокомпонентных систем, применяемом в 
плазмохимии и плазменной металлургии [1-6]. 
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(1) - уравнение, определяющее энтропию газовой 
фазы системы, где 57 - стандартная энтропия 
/-го компонента при температуре Т и давлении Р 
равном 1 атм; І щ - энтальпия /-го компонента; 
Д, - универсальная газовая постоянная; ѵ - удель- 
ный объем всей системы; щ - число атомов /-го га- 
зообразного элемента; п еі - знак и кратность иони- 
зации компонента; А,, Х е - неопределенные множи- 
тели Лагранжа; М : - содержание компонентов газо- 
вой фазы, моль/кг. 

(2) - уравнение, определяющее энтропию компо- 
нентов, находящихся в конденсированном (твер- 


дом и жидком) состоянии и образующих отдельные 
фазы, где ІІ „- внутренняя энергия конденсирован- 
ной системы. 

(3) и (4) - уравнения для определения энтропии 
раствора Х\ и Х 2 , где М щ и М„ 2 - число молей конден- 
сированных веществ, входя щ их в растворы Х\ и Х 2 , 
/V, и /Ѵ 2 - число компонентов, образующих раство- 
ры Х\ и Х 2 . 

(5) - закон сохранения массы всех химических эле- 
ментов [эА,], где Д.- число молей элементов, при- 
ходящихся на 1 кг рабочего тела; щ, п іп п м , п ]Пі - сте- 
хиометрические коэффициенты, то есть числа ато- 
мов у-го элемента в соответствующем индивиду- 
альном веществе-компоненте системы. 

(6) - уравнение состояния смеси идеальных газов, 
где Р- давление в системе. 

(7) - закон сохранения электрического заряда. 

(8) и (9) - условия нормировки для компонентов, 
вошедших в растворы X, и Х 2 . 

(10) - уравнение для определения энтропии 5 всей 
термодинамической системы. 

(11) - уравнение для определения полной энталь- 
пии /всей термодинамической системы. 

(12) - уравнение, характеризующее величину пол- 
ной внутренней энергии Ь„ равновесной системы. 

Система уравнений (1-12) является трансцен- 
дентной, поэтому для определения искомых неиз- 
вестных приходится ориентироваться не на полу- 
чение аналитического решения, а на какой-либо из 
приближенных методов вычислений, например, 
метод Ньютона. 

В МГТУ им. Н.Э. Баумана разработана прог- 
рамма "Астра 3" [2], позволяющая решать данную 
систему уравнений. Эта программа внедрена нами 
и модифицирована для расчетов равновесного сос- 
тава плазменных процессов. 

В данной работе проводились термодинамичес- 
кие расчеты равновесного состава: воздушной плаз- 
мы с добавками НЩМ - ЫБЮз, №1\'0 3 , КБЮз, 
ИЬБІОз, СзБЮз в диапазоне температур 1000.. ..6000 К 
с шагом 1000 К для давлений 0,1 и 0,001 МПа. Зада- 
вались следующие соотношения исходных веществ 
в смеси воздух: НЩМ (90 % : 10 % и 70 % : 30 %). На- 
иболее характерные результаты термодинамических 
расчетов для давления 0,1 МПа приведены на рис. 
1-4. Анализ результатов показал, что при добавке в 
воздушную плазму НЩМ в равновесном составе 
появляются продукты его термического разложе- 
ния. Так, при добавке в плазму нитрата лития в рав- 
новесном составе присутствуют - 1і, Іі 2 , ПО, ЬьО, 
Ы + , Іі 2+ , Іі 2 0 + , П N ; при добавке нитрата калия - К, 
К + , КО, К 2 , К 2 + , NN ; при добавке нитрата цезия - Сд 
С§ + , С80, С8ІМ; при добавке нитрата рубидия - КЪ, 
КЪ + , КЪО, КЬМ, К Ь>: при добавке нитрата натрия - 
N3, Ыа + , N30, № 2 . 

Увеличение концентрации вводимой в воздуш- 
ную плазму добавки НЩМ в исходной смеси с 10 % 
до 30 % приводит в расчетном равновесном составе 
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к увеличению концентрации компонентов, харак- 
терных для процессов взаимодействия воздушной 
плазмы с НЩМ, в том числе к увеличению конце- 
нтрации электронов. 



I р-0.1 МПа I 

Рис. 1. Равновесный состав системы воздухНИО 3 (70:30) 


С увеличением давления в смеси с 0,001 до 0,1 
МПа в равновесном составе отмечается увеличение 
концентраций компонентов N0, 0 2 , СО,, N0,, N 2 , 
О, О , О, , в том числе компонентов, соответствую- 
щих конкретному типу исследуемого НЩМ: нитра- 
ту лития - Ы, и, ЫО, Ы 2 0, Ы 2 + ; нитрату калия - К, 
КО, К,, К 2 + ; нитрату цезия - Сх, СкІЧ, СкО; нитрату 
рубидия - КЪ, КЪО, КЫМ, КЬ 2 ; нитрату натрия - Ыа, 
1Ча 2 , №0, №1Ч. Наряду с этим, при повышении 
давления в смеси в равновесном составе в целом от- 
мечено уменьшение концентрации электронов, Ы, 
О, С, СО, и ионов - Ій + , 1Ча + , К\ КЬ + , Ск + , для соот- 
ветствующей добавки НЩМ, обусловленное ре- 
комбинацией заряженных частиц, а также химичес- 
кими реакциями в условиях избытка окислителя. 



Рис. 2. Равновесный состав системы Аг:ИЫО г (70:30) 


Зависимость концентрации компонентов рав- 
новесного состава от температуры неоднозначна. 
Из анализа зависимости концентрации компонен- 
тов равновесного состава от температуры следует, 
что концентрация молекулярного кислорода с по- 
вышением температуры уменьшается, хотя ход 
кривой существенно зависит от типа добавки 
НЩМ и давления в смеси. С увеличением темпера- 
туры отмечается увеличение в равновесном составе 


концентраций электронов, атомарного азота, кис- 
лорода, углерода. Максимумы концентраций диок- 
сида азота соответствуют температуре 3000 К, ок- 
сидов металлов - 2500 К, №0 - 3600 К (при давле- 
нии 0,001 МПа N,0 на графиках отсутствует). Мак- 
симальное значение концентрации оксида углеро- 
да достигается при 3500 К и остается постоянным 
при дальнейшем повышении температуры при дав- 
лении 0, 1 МПа, а для 0,001 МПа заметна убыль СО 
начиная с температуры 5000 К. Максимум конце- 
нтрации СО, достигается при 7=2400 К. При повы- 
шении давления заметен сдвиг максимума конце- 
нтрации оксида азота с 2800 К до 4000 К и концент- 
рации О с 3500 К до 4000. ..5000 К (в зависимости от 
типа вводимой добавки НЩМ). При анализе кон- 
центрационных зависимостей от температуры для 
Ме и МеО характерно достижение максимума кон- 
центрации при 2300 К, за исключением Ы и ЫО. 



| р-0.1 МПа | 

Рис. 3. Равновесный состав системы воздух:С$ЫО 3 (70:30) 

Характер температурной зависимости концент- 
раций Ме + практически для всех графиков совпа- 
дает с характером зависимости для концентрации 
электронов и увеличивается с повышением темпе- 
ратуры. 



| р-0.1 МПа | 

Рис. 4. Равновесный состав системы Аг.СзЫО^ (70:30) 

Учитывая то, что НЩМ используются в качест- 
ве окислителей в ПТС, был проанализирован рав- 
новесный состав по выходу молекулярного кисло- 
рода в зависимости от типа добавки НЩМ. Из ана- 
лиза следует, что наибольшая концентрация моле- 
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кулярного кислорода соответствует вводу У КО,, и 
отмечена закономерность уменьшения концентра- 
ции молекулярного кислорода по группе щелочных 
металлов от лития к цезию; концентрация же 
электронов имеет тенденцию к увеличению от нит- 
рата лития к нитрату цезия. Закономерности в из- 
менении концентрации молекулярного кислорода 
и электронов от температуры для разных типов 
нитратов щелочных металлов сохраняются при из- 
менении исходных данных - давления и концент- 
рации НЩМ в смеси. 

С целью выяснения роли компонентов воздуха 
в процессах взаимодействия воздушной плазмы с 
НЩМ был рассчитан равновесный состав аргоно- 
вой плазмы с добавками НЩМ. Как следует из за- 
висимостей концентрации компонентов равновес- 
ного состава от температуры (рис. 2, 4), концентра- 
ции 0 2 , N 7 , N0 имеют меньшее значение по срав- 
нению с концентрациями этих компонентов в воз- 
душной плазме для аналогичных исходных данных 
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(рис. 1, 3). Концентрация электронов практически 
не изменяется, а концентрационные зависимости 
для СО, С0 2 и N,0 в диапазоне концентраций 
>10 5 моль/кг отсутствуют. Из анализа данных сле- 
дует, что молекулярный азот, кислород и оксид азо- 
та выделяются также и при термическом разложе- 
нии НЩМ. 

Термодинамический анализ равновесного сос- 
тава воздушной плазмы, взаимодействующей с 
нитратами щелочных металлов, показал, что нит- 
раты металлов щелочной группы обладают хоро- 
шей окислитель-генерирующей способностью. 
Причем, в качестве окислителей для ПТС, приме- 
няемых с целью излучения в инфракрасной облас- 
ти спектра, целесообразно использовать нитраты 
Іі, № и К, имеющих высокий выход молекулярно- 
го кислорода, а в случае использования ПТС для 
МГД-генераторов с целью повышения электропро- 
водности плазмы горения - нитратов Ск и КЪ, име- 
ющих более низкий потенциал ионизации. 
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